Differential atom interferometry beyond the standard quantum limit 
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We analyze methods to go beyond the standard quantum limit for a class of atomic interferometers, 
where the quantity of interest is the difference of phase shifts obtained by two independent atomic 
ensembles. An example is given by an atomic Sagnac interferometer, where for two ensembles 
propagating in opposite directions in the interferometer this phase difference encodes the angular 
velocity of the experimental setup. We discuss methods of squeezing separately or jointly observables 
of the two atomic ensembles, and compare in detail advantages and drawbacks of such schemes. 
In particular we show that the method of joint squeezing may improve the variance by up to a 
factor of 2. We take into account fluctuations of the number of atoms in both the preparation 
and the measurement stage, and obtain bounds on the difference of the numbers of atoms in the 
two ensembles, as well as on the detection efficiency, which have to be fulfilled in order to surpass 
the standard quantum limit. Under realistic conditions, the performance of both schemes can be 
improved significantly by reading out the phase difference via a quantum non-demolition (QND) 
measurement. Finally, we discuss a scheme using macroscopically entangled ensembles. 

PACS numbers: 42.50.St, 42.50.Ct, 95.75.Kk, 42.81.Pa 



I. INTRODUCTION 

Comparing the phase shifts obtained in two inde- 
pendent interferometric setups has several applications, 
prominent examples being the comparison of atomic 
clocks Q, and Sagnac interferometry to discriminate be- 
tween rotations and accelerations In addition, the 
comparison of gravitational forces at different points in 
space or for different atomic species allows to test pre- 
dictions ofpossible violations of Einstein's general rela- 
tivity ^ U- For such differential interferometers, the 
quantity of interest is encoded in either the difference, or 
the sum of the individual phase shifts. 

Especially for the measurement of inertial forces, atom 
interferometers promise high resolution. Here atom opti- 
cal elements like beam splitters and mirrors can be real- 
ized using Raman transitions between two atomic ground 
state levels |^. The accumulated phase difference be- 
tween the interferometer paths is encoded in the number 
difference of atoms in the exit ports labeled by differ- 
ent internal states. For example, this can be measured 
by state selective fluorescence detection. Such schemes 
have already been successfully implemented to measure 
inertial forces and the earth's gravity with a high ac- 
curacy 0, • In Sagnac atom interferometry, the goal 
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is to measure the phase shift which occurs when the 
laser setup (laboratory frame) is rotating relative to the 
frame of the freely flying atomic ensembles. This phase 
shift is given by cjf''^ = A-kA Vl m^t/h as compared to 
(flight _ 47rA r2/(Ac) for laser interferometers, where 
mat is the mass of the atoms, A is the oriented enclosed 
area of the interferometer, O, is the vector of angular 
velocity, and A is the wavelength of the light. For an 
atomic gyroscope working with *^Rb, the phase is 
10^^ times larger than the corresponding phase (^''^ht 
a light interferometer enclosing the same area and operat- 
ing at A = 10'^ nm. Hence, atom interferometers promise 
an enormously improved resolution for rotation measure- 
ments as compared to "classical" photonic devices. 

The Sagnac phase can be measured by letting two en- 
sembles of atoms pass through the interferometer from 
opposite sides. They obtain then phase shifts ±(/) due to 
the rotation of the laboratory frame, where the sign de- 
pends on the direction of propagation of the ensembles, 
and a common phase shift 9 due to effects such as an ac- 
celeration of the setup. Subtracting the phases of the two 
ensembles yields the desired phase 20, which encodes the 
rotation of the setup around an axis perpendicular to the 
plane of the interferometer. Collisions between the two 
ensembles in the interferometer can safely be neglected 
due to the low atomic densities of the ensembles. 

In the standard quantum limit for phase measure- 
ments in atomic interference experiments, the variance 
of the phase due to quantum projection noise is given by 
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(A(/))gQL = 1/-/V, where N is the number of atoms in a 
sample. By feeding the interferometer with non-classical 
states of atoms, the so-called squeezed states, this limit 
can be surpassed, with a fundamental bound given by 
(A0)2j = l/N'^^ the so-called Heisenberg limit |al9lllo|. 
Squeezed atomic states can be produced by a quantum 
non-demolition (QND) interaction of the atoms with a 
light beam lllil or by absorption of non-classical states 
of the hght Squeezed states of atomic ensembles 

might also be useful as quantum memory for states of 
light 

It has been shown that for such squeezed states atoms 
within an ensemble are entangled with each other |l4l |. 
In addition to this entanglement on a microscopic level, 
it is also possible to entangle macroscopic degrees of free- 
dom of two atomic ensembles with a similar interaction 
[lil H^ . This in principle enables teleportation of the 
macroscopic state of an ensemble. Furthermore, it has 
been shown recently that such macroscopically entangled 
ensembles can improve the efhciency of measurements of 
the components of a magnetic field Here, we try 

to exploit microscopic as well as macroscopic entangle- 
ment between two atomic ensembles in the context of 
differential interferometry, where we focus especially on 
an atomic Sagnac interferometer setup. 

In Section II we calculate the phase variance for a 
differential interferometer using non-squeezed coherent 
states in order to introduce our methodology. We show 
also how to include number fluctuations into the calcu- 
lations. These will turn out to be important later. In 
Section III we show that the phase uncertainty can be 
reduced by feeding the interferometer with individually 
squeezed ensembles. In Section IV we consider squeezing 
of a joint observable of both ensembles. In Section V we 
discuss how decoherence affects the interferometer, and 
in Section VI we compare the variances of the schemes 
discussed so far for realistic parameters. Finally, in Sec- 
tion VII the use of macroscopically entangled states in 
the interferometer is discussed. 



II. COHERENT INPUT STATES 

For a single atom we define a pseudo spin- 1/2 through 
two ground state atomic hyperfine levels |1) and |2) 0,01 
by introducing the relevant spin operators as 

= i(|l)(l|-|2)(2|) , a. = i(|l)(2| + |2)(l|) , 



= ^(|1)(2|-|2)(1|). 



(1) 



We will subsequently only consider the collective spin 
J — (*) ^'^^ ttie first, and in analogy L for the 

second ensemble; Nj^^ are the number of atoms in the 
ensembles J and L, respectively. The measurements that 
can typically be performed in atom interferometry are 



population measurements of the levels |1) and |2) by flu- 
orescence techniques 0|. Hence, only the z components 
of the collective spin vectors can be measured directly. 



II. A. Description of the interferometer 

Initially, all the atoms are assumed to be prepared in 
the state |1), leading to the following expectation value 
and variance of the collective spin vector: 

(J) = ^i, (AJ.)2 = 0, (A4,)2 = ^, (2) 

Analogous values are obtained for the ensemble L. Corre- 
sponding to the definitions in Eq. the z axis denotes 
the difference of the number of atoms in states |1) and 
1 2), whereas the phase difference between these two states 
is encoded in the x — y plane. The uncertainties stem 
from the single particle uncertainties, and correspond to 
a state with a fixed number of atoms. 

A typical interferometer sequence used for the mea- 
surement of inertial forces consists of three atom-light 
interactions as shown in Fig. ^ The first beam splitting 
Raman pulse transfers all the atoms from the ground 
state |1) to the superposition -^(jl) + 12)). Atoms trans- 
fered to state |2) obtain a momentum kick of two photon 
recoil if the two Raman lasers are counter-propagating, 
so that the partial waves delocalize, as depicted in Fig.^ 
The second pulse exchanges the populations, |1) ^ |2), 
and deflects the partial waves. Finally, they are recom- 
bined in the last interaction zone, acting as a beam split- 
ter. These pulses can be represented as rotations of the 
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FIG. 1: Scheme of the atom interferometer. ^ and vr label 
the beamsplitting and the mirror pulses, respectively. The 
populations in the two exit ports are detected state-selectively 
via fluorescence measurements. 



collective spin vectors (J) and (L) around an axis in the 
x — y plane, the angle being given by 7r/2 for beam split- 
ters and by tt for mirrors. For a fixed coordinate system 
the angle between the x axis and the rotation axis is given 
by the laser phase The laser phases change if the 
setup (laboratory frame) rotates with respect to the path 
of the freely flying ensembles, which causes the Sagnac 
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phase shift. This change corresponds to a rotation of the 
collective spin vectors in the x — y plane around the z 
direction. 

In order to make the scheme applicable to a more gen- 
eral scenario for differential interferometers, we model it 
as follows for each of the ensembles labeled by J and 
L, respectively: the first beam splitter rotates each col- 
lective spin vector by 7r/2 around the y axis, then we 
collect all the phase shifts occurring in the interferom- 
eter in rotations around z of J and L by $j and <I>l, 
respectively. Finally, a tt and a 7r/2 pulse, both around 
X, implement the mirror and the final beam splitter, re- 
spectively. These last two pulses can be combined to a 
single rotation around the x axis by — 7r/2. All relevant 




FIG. 2: (Color online) Interferometric scheme for a coherent 
input state: (a) collective spin after the first beam splitter, 
(b) the total phase 4> accumulated in the interferometer results 
in a rotation around z, before (c) a final rotation around x 
by — 7r/2 implements the mirror and the final beam splitter, 
encoding the phase in the z direction. 



L, because the Sagnac phase takes a different sign, 
depending on direction in which the ensemble passes the 
interferometer. We define the phase operator for the case 
of coherent states (cs) as 



Nj Nl 
and obtain to the first order in ch and 6 



(A0C 



1 



1 



(5) 



(6) 
(7) 



The variance of (j)cs corresponds to the standard quan- 
tum limit. Note that 9 can be obtained in an analogous 
way. From the definition it is obvious that determining 
the number of atoms is important for the calculation of 
the phase shift. A major source of error will come from 
fluctuations in the number of atoms of the ensembles and 
hence we will discuss how to include this process into the 
calculations in the next section. 



II. B. Number fluctuations 



interferometric steps are described in Fig. |3 where the 
atomic spin vector is depicted together with a disc repre- 
senting the corresponding uncertainties of the spin opera- 
tors. To simplify notation we will take the state after the 
first beam splitter as the initial state: J'" = 7?.j,(7r/2)J, 
see Fig. EJa). Here TZi{a), i S {x,y,z}, is the matrix 
rotating a vector by the angle a around the direction e^. 
In the Heisenberg picture, and neglecting collisions be- 
tween the atoms within the ensemble, the spin operator 
changes according to J°"* = 7^:^(-7^/2)7^4$J) j'". This 
leads to 



J^™ cos $ J - J™ sin 

7in 

•J z 

-Jj"sin$j-Jj" cos$j 



(3) 



It is also possible to consider a balanced atom interfer- 
ometer in which all rotations are around the x axis. In 
this case, an extra 7r/2 shift around z leads to the same 
result. The advantage of the extra pulse is that for the 
initial state of Eq. ||2J) 

(ir) = -^Bin$,«-^$, (4) 

for small angles $j, while with the unmodified balanced 
scheme ( J°"*) would be proportional to cos $j, and hence 
sensitive to $j only in second order around <i>j = 0. 

As explained in the introduction, the phase will be 
given by $j = <\) + for J and by $i = —0 -f Q for 



There are two sources of deviations of the number of 
atoms: the preparation process and the number measure- 
ment process. The atomic ensembles produced from the 
source are best described as a statistical mixture of states 
with different atom numbers, but we will assume that the 
final number measurement projects onto a number state 
with N J and Nl atoms in the two ensembles. Defining 
iV = {Nj + Nl)/2, we assume that \Nj - Nl\ = jVn, 
which reflects the variances of the number operators Nj^l 
of the ensembles after the production. Thus, 7 is the pa- 
rameter which describes how well the atom numbers in 
the two ensembles match. 

We treat the number measurements by introducing op- 
erators SN with 



{SN) = 0, [A{SN)f = aN, 



(8) 



where a describes the quality of the number measure- 
ment. For fluorescence measurements, is given by 
the mean number of times an atom goes through the flu- 
orescence cycle and scatters a photon which subsequently 
is registered in the detectors |20j . Typical values in nowa- 
days experiments are « 50 . . . 100. We replace Nj by 
N^+SN^P+SN^'' and J™* by J°''''° + {SnP -5Nf'')/2. 
Nj refers to the actual number that would have been 
measured in perfect number projection measurements, 
and in SNj^"^ the index i corresponds to the atomic levels 
\^). 
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With these substitutions, we obtain 



(9) 



where terms of higher order have been neglected. The 
contribution from the number fluctuations is in agree- 
ment with Ref. [i^- The expectation value of (^cs is 
shifted for a ^ 0. As this shift is of the order of the sec- 
ond term in (A(/)cs)^, it is of second order of the standard 
deviation only, and can thus safely be neglected. From 
the expressions it is clear that a <C 1 is required in order 
to reach the fundamental limit for the phase resolution 
using coherent input states. 

In the following, we will drop the superscript on the 
atom number, as long as there is no danger of confusion. 



III. SEPARATELY SQUEEZED ENSEMBLES 

It is known that by taking squeezed input states it is 
possible to surpass the standard quantum limit in inter- 
ferometry (Qj] . In this section, we consider the case where 
both ensembles are squeezed separately with the method 
introduced in i.e., by a QND interaction with a laser 
beam shortly after the first beam splitter. 



III. A. Squeezing a single ensemble 

We assume to have the situation of Fig. O i.e., the 
electromagnetic field mode ai couples states |1) and |3), 
and a2 couples states |2) and |4). Here transitions |1) ^ 
1 4) and |2) ^ |3) have to be suppressed to the first order 
in the coupling constant f|lll|. see also the discussion in 
Section V). For the light, an effective spin vector can be 



|3)- 



|4>- 



|2>- 



FIG. 3: For the QND interaction, the levels |1) and |2) are 
coupled ofF-resonantly to states |3) and |4). 

defined, the so-called Stokes vector. Its components are 
given by 



5. 



- (^a\ai - ala2^ , Sx ^ ^ (^a|a2 -I- alai^ ,(11) 

Sy = ^(a|a2-4ai), (12) 

where a\ and create a photon in mode ai and a2, 
respectively. Sz measures the difference of photons in the 
two modes. By an appropriate choice of the parameters. 



the Hamiltonian describing the interaction of the light 
with the first atomic ensemble can be brought to the 
form [ni 



H = hQjSzJz 



(13) 



with a frequency flj. Due to the interaction, the collec- 
tive spin vector of the atoms is rotated around the z axis 
by xj^z, with the atom-photon coupling XJ — s-nd 
the effective interaction time t. The Stokes vector un- 
dergoes the same variation with XJ^z replaced by xjJz- 
This rotation of the Stokes vector is due to the Faraday 
effect of the light passing the atoms, while the rotation of 
the atomic spin vector is due to an AC Stark shift origi- 
nating from the light field. The coupling x/ is given by 

E3 



XJ 



ciF2 



A 



A2' 



(14) 



where A and L are cross section and length of the atomic 
sample, F and d are line-width and dipole moment of the 
atomic transition, respectively, and A is the detuning 
from the atomic resonance frequency w. 
After the interaction. 



5™* = sin{xjJz)S'^ + cos(xj J.)4' 



and if initially 



yx, {AST? 



0, (A^^ 



(15) 



where nj is the number of photons in the ensemble J, 
then we can effectively replace 5^'" by its macroscopic 
expectation value nj/2. Furthermore, developing the 
trigonometric expressions and assuming that NjXj 1, 
we obtain in leading order 



S° 



njX 



J" 



SI 



V ■ 



(17) 



Hence a measurement of the y component of the outgoing 
light vector gives information about {Jz), while Jz itself 
is not affected by the rotation around the z axis. 

If such a QND measurement is performed after the first 
beam splitter of the interferometer, then the operator 



j ^ — j ^ 



njXJ 



-S?, 



(18) 



measures the difference between the z component of the 
ensemble's atomic spin vector after the second beam 
splitter and the estimated value after the first one. The 
fiuctuations of this operator are reduced as compared to 
J°"* [nj, while the fluctuations of Jy are enlarged, which 
is depicted in Fig. 21 Hence, the state of the atomic spin 
vector is squeezed in the z direction. 
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III.B. Modified interferometric sclieme 

We modify the scheme introduced in Section II by in- 
serting the QND interaction shortly after the first beam 
sphtter, foUowed by an extra rotation around x by 7r/2, 
which rotates the uncertainty elhpse such that the phase 
uncertainty is reduced as desired, cf. Fig. 01 In the 
experiment, the latter pulse must not transfer momen- 
tum to the particles, which can be achieved by using 
co-propagating Raman lasers for this step, provided that 
the two transition frequencies are approximately equal, 
so that the two recoil momenta cancel in the transitions. 



provided that a is small enough, the variance is dom- 




FIG. 4: (Color online) The interferometric scheme for 
squeezed input states: (a) collective spin after the first beam 
splitter, (b) the QND measurement prepares a spin squeezed 
state which (c) is rotated around x by 7r/2; step (d)^(e) as 
step (b)-+(c) in Fig.H 

The outgoing spin vector now is calculated as J°"* — 
7^,(~7^/2)7^^($J)7^,(^/2)7^,(x^,)J"^ and in analogy 
for the second ensemble L, with corresponding Stokes 
vector T. 

In comparison to Eq. ||SJ|, J°"* now has to be corrected 
to incorporate the spin squeezing (ss) as described above: 



— _ ^ ( TOUt 



njXJ 
2 

nhXL 



-SI 



(19) 



Calculating the expectation value and variance as before 
yields 



nx 



4- 



1 

Nl 



+ a( — 
\Nj 



NlJ 



(20) 
(21) 



Here it has been assumed nj = =: n and XJ — XL ='■ 
X- Deviations from these assumptions enter the variance 
only in higher order terms as long as ||xj| — |xl||/(Ixj| + 
I XL I) ^ 1, and similar for nj. Further assumptions 
leading to these expressions are ^x^ <C 1 as mentioned 
before, as weU as nx^ < VS> ■ [\/iV(6'^ -I- 0^)]"^ Now, 



inated by the first two terms scaling as N 



J/L- 



They 



originate from the projection noise of the light, and in 
principle allow to improve the resolution below the stan- 
dard quantum limit. However, nx^ equals, except for 
a factor of order unity, the fraction of atoms which are 
lost due to spontaneous processes during the squeezing 
process, cf. Section V. Thus, nx^ <C 1 is necessary, 
and even though we obtained a Heisenberg-like scaling 
(A(/)ss)^ ~ l/Nj we are far from reaching the Heisen- 
berg limit. 

Furthermore, as it becomes clear from the second term, 
the resolution is limited by the accuracy of the fluores- 
cence number measurements. These measurements are 
necessary in any case to determine the phase, because 
( J°"*) and (-fc™') have to be rescaled properly. However, 
Jz and Lz itself can be measured using another QND 
interaction. As will be shown in the next section, this 
reduces the dependence on the quality of the number 
measurements. 



III.C. QND output measurement 

Let us consider now a modification of the scheme us- 
ing squeezed states, where a second QND laser beam is 
sent through each ensemble shortly after the last beam 
splitter. We define 



1 2 

1 2 

Nl nx 



f QOUt oout \ 



_| / rpOUt 

-» r \ y.r 



rpOUt 

y 



(22) 



where the extra index -I- is supposed to indicate the ad- 
ditional QND measurement. Furthermore, and 
correspond to individually prepared light pulses used to 
read out Jz and L^, respectively. The resulting expecta- 
tion value and variance in leading order are 



at 



nx 



1 1 
2/11 

2 



Nj 



Nl 



1 



mj ANl 



(23) 



(24) 



Let us compare the modified scheme including an 
additional QND measurement with the scheme using 
squeezed states. The magnitude of a in the variance of 
the modified scheme is effectively reduced for small 9 and 
0, a — > a{9'^ + </>^)/4, as compared to Eq. (|^ . Hence, 
the dependence on the number measurement is reduced. 
Furthermore, the leading term shifting the expectation 
value from the desired result (f> is smaller by a factor 9 as 
compared to Eq. H2U|) . 
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As a further possible advantage, the effect of a non- 
symmetric atom-hght interaction is compensated if the 
couphng of the read-out Q ND p ulse is similar to the cou- 
pling of the squeezing pulse j2l| . However, this advantage 
is probably not very relevant for the Sagnac interferom- 
eter considered here due to the long time-of-flight of the 
ensembles in the interferometer between the two QND 
pulses. 

A disadvantage of this scheme is that the first term of 
the variance of Eq. H24() comes with a factor of 2 because 
of the two projection measurements of the light necessary 
per atomic ensemble. It is possible, although technically 
demanding, to reduce this contribution by re-using the 
light from the first QND interaction for the read-out. In 
this case, x ^ is needed in the second interaction 
in order to obtain the difference J™* — J™ |2l|, and in 
analogy for the second ensemble. The sign can also be 
achieved by a tt rotation of the atom spin vector around 
the X axis in between the final beam splitter and the 
QND read-out pulse. 



We define the phase operator for the scheme employing 
jointly squeezed (js) ensembles as 



1 



1 



-L° 



^ g out 

nxN ^ 



(27) 



Note that in this definition and L°"* are divided 

by the atom numbers of the respective ensembles as be- 
fore, while the QND measurement yields an estimate of 
J^°"' -I- L™' without such a correction, cf. Eq. As 
a consequence we expect that we loose the advantages of 
squeezing if the numbers of atoms in the two ensembles 
differ strongly, i.e., if 7 3> 1. We find in leading order 



nx 



2a 
1 2a 

w 



27V2 



(28) 
(29) 



The importance of the number fluctuations at the prepa- 
ration stage is reflected in the fact that in order to arrive 
at these equations, the assumption 



IV. JOINTLY SQUEEZED ENSEMBLES 



7V^(nx ) +r) < 1 



(30) 



In the preceding section we have seen that the 
1 / {nx'^ Nj) term in the variance comes with a factor given 
by the number of QND interactions with different ensem- 
bles of light, cf. Eq. H24I) . For this reason let us consider 
the case of preparing the initial state of the two atomic 
ensembles with only a single QND pulse that interacts 
with both ensembles consecutively. The first interaction 
(with ensemble J) transforms j'" 7ez(xS'f ) J'", the 
second interaction (with ensemble L) transforms L'" 
7?.z(x5'l,")L'", because S™ itself remains unchanged dur- 
ing the QND interaction. The Stokes vector S transforms 
as S°"' = 7^^(x^i")7^^(xJf )S'", and the y component of 
the outgoing light is given by 



cos(x4") cos(x )^^" + sinixJhS^ 



such that for iVx^ 
Eq. Hlt)|) . we have 

^out ^ 



cosixJnS^ - sin( Ji"x)^r (25) 



<C 1, and S initially prepared as in 



nx 



L" 



(26) 



Measuring 5°"* thus reveals information about J™ + 
and performs a squeezing operation on this joint opera- 
tor. Now we apply the same operations as before to the 
ensemble J, but for the ensemble L we perform a rota- 
tion by TT around the x axis before the final measurement. 
After the extra pulse, the Sagnac shift cf) is effectively en- 
coded in the sum of the z components instead of into the 
difference. 



is necessary in addition to the assumptions leading to 
Eq. 1)21(1 . Furthermore, now 1/ {nx'^N'^) is the leading 
term only if nx^7^/8 ^ 1. 

A QND measurement could also be used after the in- 
terferometer to directly read out the joint observable 
J°^* + Lf* by defining 



2 

nx N 



-(s° 



s° 



(31) 



where again the indices r refers to the read-out QND 
measurement. Notice that in this way it is not 
possible to measure the correct rescaled observable 
J°^*/iVj -I- L°^^ /Nl, and consequently there is an impor- 
tant contribution from the difference of the atom num- 
bers in the expectation value already: 



Nl-Nj ^ a 
Nl + Nj 2N 



^. (32) 



Compared to the variance of the scheme using jointly 
squeezed ensembles without the QND read-out measure- 
ment, Eq. (|29|l . the dependence of the variance on the 
number measurements and on the atom number differ- 
ence in the two ensembles is reduced for small 6, 6: 



(A0, 



nx^N"^ 

v,2 



a 

2N 



7_ 
'8iV2' 



(33) 



However, the 7-dependent correction to the expectation 
value {(t>is+) is only negligible compared to the standard 
deviation A0js+ if -f'^Nnx'^0'^ /8 <C 1. This is generally a 
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stronger criterion than the limit on 7 encountered with- 
out the QND read-out, cf.Eq. (j^ . 

The offset can be compensated by using an estimate 
for 6 from the final fluorescence measurement 



_ TOUt 



1 

'n'l 



-l: 



to define a corrected phase operator 



Nl + Nj 



which takes into account the bias of 1 
to leading order 



(34) 



(35) 



Wc find that 



acj) 
^ 2N 
2 a 



2m 

2 



nx^N'^ 2N 8iV2 



(36) 
(37) 



Hence, the 7-dependent bias in the expectation value is 
reduced by a factor a^/ {N^/'^d) <^ 1, while the 1/N term 
in the variance still has a factor acjP. 

These are the same advantages that we also found 
in the case of separately squeezed ensembles with QND 
read-out, cf. Eq. H24|l . However, the contribution of the 
term proportional to l/nx^N'^ is reduced by a factor of 2 
in (A(^?j,^)^, because only two instead of four projective 
measurements of the light are necessary in this case. 



V. DECOHERENCE 

The attainable squeezing is limited by the absorp- 
tion of photons during the interaction between light and 
atoms ^22] . Each atom which absorbs and subsequently 
spontaneously emits a photon is no longer correlated to 
the rest of the atoms, but still adds to the variance. We 
estimate the number of atoms contributing to such an 
uncorrelated background as the number of scattered pho- 
tons nn, where k — A^j^F/A is the optical density. Then, 
in the limit of k <C 1 and for just a single ensemble, one 
finds 22] 
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for the collective atomic spin vector and 
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for the Stokes vector. The leading order corrections to 
the variance (A(/)ss)^ given in Eq. (|21|l for the case of 
separately squeezed ensembles reads 
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and similarly for the other schemes. We will use this 
estimate in the following discussion and leave an in-depth 
analysis of decoherence processes, e.g., following the lines 
of 22], to further investigations. Generally, according 
to Eq. (|38l) also the expectation value changes due to 
decoherence. This can be accounted for by rescaling, and 
doing so gives only higher-order corrections to Eq. (|42|l . 

For usual choices of parameters, the last term in Eq. 
()42|l is negligible, while the contribution proportional to 
is comparable in size to the terms in Eq. H21|l . This 
limits the coupling x E^nd thus the achievable squeezing. 
For all other experimental parameters fixed, there exists 
an optimal choice for the detuning A and thus for x [see 
Eq. (|14|) ] which minimizes the variance. Taking into 
account only the first term in Eq. I|21|l . and in the limit 
A ^ F, this optimal choice of A leads to a minimal value 
for (A0,,)2: 
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The scaling is thus no longer as in the Heisenberg limit, 
but it is still better than in shot-noise limited measure- 
ments (see also 113] )■ 

In addition, the decay of the states |1) and |2) dur- 
ing the interferometer step has to be taken into account. 
Choosing long-lived hyperfine ground-state levels to im- 
plement |1) and |2) minimizes this decay. Also, spin 
squeezed states have been shown to be robust with re- 
spect to both, particle loss and dephasing in con- 
trast to, e.g., GHZ states, which are maximally fragile 
under particle losses. 



VI. COMPARISON OF THE SCHEMES 

To analyze the performance of the schemes discussed 
in the preceding sections we will fix the number of pho- 
tons as n = 10^^ and take N — 10^" as a reasonable 
parameter for the mean atom number per ensemble. We 
will first consider a close-to-ideal scenario and assume 
that the noise from the fluorescence measurements can 
be neglected by setting a = 2 x 10~^. Also, we will 
set 7 = 10, which for N = lO^*' atoms corresponds to 
\Nj — Nl\ = lO^^TV, and we will initially not include 
decoherence. Fig. O (a) shows the scaling with N of 
(A(/))2/(A(/)cs)^, i.e., of the various variances normalized 
to the case of coherent ensembles. For all the methods 
involving squeezing of some observable a Heisenberg like 
scaling is visible. The offsets of these curves are given 
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FIG. 5: Double-logarithmic plot of the variances {A(j))^ for the methods discussed in the text normalized to the variance for 
coherent states (Ai/ics)^, as a function of the mean number of atoms N. The scale on the right hand site gives the noise reduction 
in dB. In both figures n — 10^^, and we fixed 9 = (j> — 0.01 in order to operate close to the point of maximal sensitivity of 
the interferometer. In (a) we consider the close-to-ideal scenario with a — 2 x 10~^, 7 = 10, and decoherence is not included. 
X ~ 3.23 X 10^^" is used, corresponding to a detuning A = 2.28 x 10^" s~^ for the Rb D2 line and a cross section A — 0.3 mm^ 
of the laser beam. In this close-to-ideal scenario, the scaling as iV~^ is visible for all the methods employing squeezing. The 
graphs for separately squeezed ensembles (without QND read-out) and for jointly squeezed ensembles with QND read-out and 
corrected expectation value lie on top of each other. In the case of a joint QND read-out, failing to correct the expectation 
value for the contribution from 6 results in a scaling as 1/N for iV > 10^. (b) Realistic scenario with a = 2 X 10"^, 7 = 10*, 
and including decoherence. In each case and for each value of N the detuning has been adjusted to minimize the variance 
[2^. The inset shows the corresponding optimal values of the coupling parameter x- For large atom numbers it is clearly seen 
that the schemes which do not employ a final QND measurement are strongly afi'ected by the limitations from the fiuorescence 
detection. The curve for jointly squeezed ensembles using a joint QND readout lies outside the range of the figure. 



by the numbers of QND measurements performed, i.e., 
by the factor multiplying the l/{nx^N'^) term in the 
variance of each of the considered schemes that involve 
squeezing. For the measurement of (j) via reading out 
Jz + through a QND interaction, the term propor- 
tional to 6 shifting the expectation value [see Eq. (|32|l ] 
has been included into the variance, in order to allow for a 
fair comparison. It is this term which makes the variance 
scale only proportional to 1/iV for N > 10^. Obviously, 
correcting this contribution of 9 as described in Section 
IV avoids this term and maintains the improvement by a 
factor of 2 compared to the case of squeezing and QND 
measuring both ensembles separately. 

In Fig. |S1 (b) , the relative variances are plotted for re- 
alistic experimental parameters and including decoher- 
ence. a = 2 X 10~^ corresponds to 50 fluorescence cycles 
per atom, and 7 = 10* is equivalent to a difference of 
the number of atoms in the two ensembles of 10% of the 
mean number iV at iV = 10^°. With all other parame- 
ters fixed, for each value of N the interaction strength x 
is determined by choosing the detuning A from atomic 
resonance such that the variance (A0)^ is minimized, see 
Eq. gSl) ^25,]. 

As can be seen from Fig.|Sl(b), in such a realistic sce- 
nario the noise reduction obtained from squeezing de- 
creases for all methods. For large atom numbers, in all 



cases the variance scales as 1 /N due to decoherence and 
the noise from the fluorescence measurements. For all the 
procedures not involving a QND read-out, the strong in- 
fluence of the latter contribution can be observed, though 
for N = 10^° atoms the total noise still is reduced by 
around 7 dB with respect to the limit set by quantum 
projection noise. On the other hand, the read-out via 
QND measurements allows to reduce the noise by more 
than 10 dB compared to this limit, while this method 
does not require an additional experimental setup as 
compared to the scheme where QND measurements are 
performed only on the incoming atomic ensembles. 



While in the closc-to-ideal case squeezing of a joint ob- 
servable gives an advantage of a factor of 2 in the variance 
compared to individual squeezing and read-out (corre- 
sponding to a 3 dB noise reduction), for an experimen- 
tal reasonable scenario this advantage reduces to about 
1.5 dB. Of course the method of squeezing a joint ob- 
servable of both ensembles has a basic advantage since 
it only needs a single squeezing operation instead of two, 
and thus less technical effort is necessary. 
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VII. ENTANGLED ENSEMBLES 



These measurements yield both angles: 



Julsgaard et al. demonstrated experimentally in 
Ref. [lal the generation of macroscopic entanglement be- 
tween two atomic ensembles. The scheme to generate 
such a macroscopically entangled state, described first in 
Ref. is motivated by the fact that under the ideal 
condition of 7 = two commuting joint observables can 
be constructed from J and L: 

{[Jy^LyJ,+L,]) ^{Nj-Nl)^0, (44) 



i.e., Jy — Ly can be measured without affecting 
and vice versa. This can be seen directly from 



(7^,(x5,)J'"-7^,(x5,)L" 



\ _ Tin 
)y — 'Jy 



L" 



(45) 



i.e., the first QND interaction leaves the difference of 
the y components unaffected. Thus after squeezing the 
sum Jz + Lz, also the difference Jy — Ly can be squeezed 
without loosing the information gained in the first mea- 
surement. To realize this experimentally in the interfer- 
ometer, after the first squeezing interaction J is rotated 
by a classical 7r/2 pulse around the x axis so that Jy — > Jz 
while L is rotated by — 7r/2 around x giving —Ly — s- L^. 
Then a second laser pulse, prepared again as in Eq. (|f 6|l . 
interacts consecutively with both ensembles and thus fi- 
nally carries information about Jy ~ Ly. The outgoing 
state corresponds to a macroscopically entangled EPR 
state It is now a natural question to ask whether 

entangled atomic ensembles are of use in Sagnac atom 
interferometry. 

For the schemes discussed so far, the collective spin 
vectors lie in the x — z plane after the last step of the 
interferometer. Therefore always {Jy — Ly) =0, and an 
additional operation is necessary in order to encode phase 
information in the y components as well. This can be 
achieved by rotating both ensemble vectors J and L by 
an angle if and —ip around the x axis before and after 
the interferometric phase is applied, respectively. In this 
way the plane of rotation of the phase shift is effectively 
tilted by ip around the x axis. 

The measurement process now consists of first rotating 
J and L by ±7r/2 and using another QND interaction to 
measure the sum of the z components. This measure- 
ment, scaled correctly by a iy9-dependent factor, reveals 
(j). To be more explicit, the corresponding operator for 
entangled ensembles (EE) reads 



1 



See 
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( ^™ 
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(46) 



A measurement of the sum of the y components can be 
realized after another rotation around x by either a QND 
or a projection measurement. In the former case 
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where the offsets can be corrected as above. For parame- 
ters as in Section V. the leading terms of the correspond- 
ing variances read (at 6* = and (f> = 0) 



(A<^ee)' 
(A^ee)' 
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cos^ ip N^nx^ 
1 2 



' ip N^nx^ 
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(50) 
(51) 



Changing ip allows to trade in a lower variance of one 
component for a higher variance of the other. But these 
variances are only scaling with if 7 is close to zero, 

otherwise the last term cx j'^nx^ /N in Eqs. H50|l and H51() 
is dominating the scaling. However, 7 « is an obvious 
requirement in this case, as otherwise the commutator 
does not vanish in Eq. H44|) . and thus the two squeezing 
operations are not compatible. 



VIII. CONCLUSION AND OUTLOOK 

We have presented and compared in detail several 
methods to improve the detection of a differential phase 
shift of two atomic interferometers beyond the standard 
quantum limit, having in mind especially the application 
to Sagnac interferometry. For this purpose, we have an- 
alyzed the squeezing of individual and joint observables 
and, in both cases, the read-out of the interferometer via 
fluorescence detection of the atoms only or by an addi- 
tional QND interaction. 

If decoherence and measurement imperfections are ne- 
glected, all the methods of squeezing improve the be- 
havior of the variance of the differential phase to a 1/ N'^ 
scaling, modified by a factor k/{nx^) ^ 1, which is deter- 
mined by the number k of QND interactions involved, by 
the number of photons n, and by the coupling x between 
atoms and photons. In the case of jointly squeezed ob- 
servables, we found that this limit can only be attained 
if some constraints on the difference of the number of 
atoms in both ensembles can be fulfilled. In all cases, 
the achievable squeezing is limited by decoherence due 
to the absorption of photons during the QND measure- 
ment. 

Using fluorescence measurements to read out the 
atomic spins after the interferometer always produces 
additional noise scaling as 1/N due to the photon shot 
noise. As an alternative method, a QND measurement 
can be employed to read out the final state of the in- 
terferometer. Although in this case fluorescence mea- 
surements are still necessary to determine the number 
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of atoms in the two ensembles, their contribution to the 
noise is reduced to a large extent. We have shown that 
the best method to achieve this is to perform squeezing 
and read-out via a QND measurement of a joint observ- 
able of the two ensembles, provided that the difference 
between the number of atoms in the two ensembles can 
be made smaller than approximately 10% of the mean 
number of atoms. This procedure minimizes the number 
of QND interactions necessary, thereby minimizing the 
factor multiplying the l/N"^ term in the variance, and it 
reduces the experimental effort. 

Finally, we considered the creation of a macroscopi- 
cally entangled state of the two atomic ensembles via 
squeezing of two non-local, commuting observables. We 
showed that in this case both the sum and the differ- 
ence of the phase shifts can be measured with a variance 
scaling with and that the relative uncertainty can 

be shifted between both quantities. However, this scaling 
can only be reached here if the number difference between 
the two ensembles can be made very small. Therefore, 
it would be desirable to identify methods to control the 
numbers of atoms in the ensembles, for instance by em- 
ploying the superfluid - Mott insulator transition in an 
optical lattice embedded in a weakly confining harmonic 
potential jl^ H^]. Controlling this confinement, which 
plays the role of a local chemical potential, the number 
of atoms in the Mott phase at T > can be controlled 
and should only depend mildly on the total number of 
atoms in the system. A detailed analysis of this idea is 
left for future investigations. 
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